Low energy excitations in fermionic spin glasses: 
a quantum— dynamical image of Parisi symmetry breaking 
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ON : 

We report large effects ol Parisi replica permutation symmetry breaking (RPSB) on elementary 
excitations of fermionic systems with frustrated magnetic interactions. The electronic density of 
| states is obtained exactly in the zero temperature limit for (K — l)-step RPSB together with exact 

relations for arbitrary breaking K, which lead to a new fermionic and dynamical Parisi solution at 
K = oo. The Ward identity for charge conservation indicates RPSB-effects on the conductivity in 
£ — , metallic quantum spin glasses. This implies that RPSB is essential for any fermionic system showing 

spin glass sections within its phase diagram. An astonishing similarity with a neural network problem 
1 , 1 is also observed. 

« ! 

PACS numbers: 71.20.-b, 71.55.Jv, 75.10.Nr 

We present the first solution to the question whether and in which way Parisi replica permutation symmetry break- 
ing (RPSB) and the related nonconstant part of the Parisi spin glass order parameter function q(x) |^||] are displayed 
in the low temperature many body theory of fermionic systems with frustrated Ising-interactions, emphasizing the 
T = 0-limit in particular. The Parisi function q(x), defined on the interval < x < 1 is known as the apparently 
exact solution of the infinite range classical spin glass. The x-dependence is known to be comparable with a non- 
\ trivial time-dependence, induced by Glauber dynamics Q, of the spin autocorrelation function < o"(t x )<j(0) > (small 
x corresponding to large times). As shown by Parisi the function q{x) assumes a plateau value within x\ < x < 1 
and differs from this sofar conventional single order parameter picture only within < x < xi = 0(T), where it 
, decreases towards zero at x — in the absence of a magnetic field. The RPSB-effect appears to disappear with 

T — > 0. Nevertheless we find and report here a large 0(T )— effect to persist in many important physical quantities 
t-H , of the fermionic Ising spin glass, which is a minimal quantum generalization of the classical Sherrington Kirkpatrick 
model. This includes replica-diagonal fermion Green's function and fermion density of states, where at any step K 
of RPSB the set of different order parameters is seen to determine the quantum-dynamical behaviour of the fermion 
propagator and of vertex functions. These effects are complementary to and not in contradiction with recent replica- 
symmetric descriptions of T — quantum spin glass transitions Q. Parisi-RPSB jjj is seen to decide the qualitative 
and quantitative features of the low energy excitation spectrum. While results are presented for an insulating model, 
the effect appears to be rather model-independent and should hence be felt in transport properties of models with 
&\ , additional hopping hamiltonian for example. 

' The presence of spin glass phases within phase diagrams of interacting many-fermion systems such as 
HighT c superconductors, heavy fermion systems, and semiconductors are nowadays recognized with increasing at- 
tention ^-0. Many of their characteristic properties cannot be answered by considering these phases as isolated 
magnetic phenomena, which means that their common origin, their coexistence and competition with charge-related 
phenomena, and even far-reaching links into other fields of physics must be understood in terms of fermionic rather 
than spin space models. The relationship between conductivity behaviour and magnetism, mainly antiferromagnetism 
, up to now, has acquired a prominent place in the conscience of theorists and experimentalists, due to the remarkable 
progress in the field of strongly correlated fermion systems during recent years JsJJ^] . 

In this Letter we wish to provide results which evidence the fact that fermionic spin glasses also link closely glassy 
magnetic order and transport behaviour; further similarities between Hubbard model and the fermionic spin glass 
have been traced back to the particular role of the Onsager-Brout-Thomas reaction field m for all these systems, 
as can be observed by comparing Hubbard-CPA- || with fermionic TAP-equations . 

Spin- and charge-excitation spectra of fermionic spin glasses must be evaluated in order to construct a meaningful 
many body theory. This Letter focusses on the effect of Parisi replica permutation symmetry breaking (RPSB) on the 
single fermion density of states (DoS), hence on the fermionic Green's function, and, by virtue of the Ward identity 
for charge conservation, also on vertex functions, thus on the entire ensemble of quantities that provide the basis of 
many body theories for fermionic systems with frustrated interactions. 

It is known that replica-diagonal quantities like the linear equilibrium susceptibility x feel Parisi symmetry breaking 
even at T = despite the fact that the nontrivial part of the Parisi function only lives on an interval of width 
T. 

The susceptibility had been analysed by Parisi for the standard SK-model. He found a rapid convergence towards 
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the exact result as the number of order parameters increased, this number being equal to K + 1 in the SK- and 
equal to K + 2 in fermionic models. While the low temperature regime of the SK-model had not been of particular 
interest from the point of view of phase transition theory, it becomes highly important for fermionic spin glasses, 
since the T = 0-theory of excitation spectra plays a crucial role and, for the additional reason that some models 
exhibit quantum phase transitions along the T = 0-axis. Parisi nevertheless analysed the low T regime @] of the 
classical SK-model finding that K-step RPSB on one hand provided increasingly good approximations but failed to 
completely remove the negative entropy and the instability problem at low enough temperatures unless K — > oo. 
In this Letter the effect of one step RPSB (K = 1) on the density of states is presented in detail, followed then by an 
analytical relation valid for all K, which allows to determine the type of excitation spectrum present in the full Parisi 
solution for the fermionic Ising spin glass. Despite the fact that the regime of deviation from a replica-symmetric 
spin glass order parameter is only of 0(T), we find that it has a large O(T ) -effect on the fermion density of states, 
the one-particle-, and many-particle Greens functions at T = 0. This density of states is derived as usual from the 
imaginary-time (disorder-averaged) fermion Green's function [— < T T [a(r)a' (0)] >] a v, which is one of the decisive 
quantum-dynamical elements of any many body theory of fermionic spin glasses. This illustrates that, unlike the 
usual picture of a Parisi solution being just a static order parameter function, the fermionic picture must include the 
qualitative extension to dynamical quantities. Those become drastically altered by the nontrivial part of the Parisi 
solution which is otherwise invisible at T — 0, hence providing a quantum-dynamical image of RPSB. We first consider 
a generalized Parisi solution of the infinite-range fermionic Ising spin glass model. Its grand canonical Hamiltonian 

H = ~y] Jijdi&j - H &i - + hii), °t = n A _ %) ™w = a L a ^' (!) 

ij i i 

with fermion operators a, a' and represents a Fock space extension of the SK-model. The magnetic couplings Jij of 
this insulating model are independent gaussian variables with zero mean value. The chemical potential controls the 
occupation of magnetic and nonmagnetic states, where the latter ones reduce the freezing temperature, and leads to 
remarkable effects in the tricritical phase diagram |pT[ . We restrict our discussion to half- filling. The fermion Green's 
function can be derived as Q — j^j-lnB from the generating functional (generalization from K = 1 step RPSB, given 
here for the sake of simplicity, to arbitrary K is standard) 

r oo , (« 7 ) , '"A 2 r ™/ m aim 
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Q-i — L a-(ai-l)m+l 



with a bare propagator Go{ti) — (*£/ + m) 1 and a magnetic field H included in the effective field H({z\ a ^}) 

H + JJ2~/ — 97+i z t* 7 \ where q a = q, q K+1 — 0. Fermionic fields are denoted ip,%jj, and n,fj. Spin (decoupling 
fields z 7 , carrying_ a Parisi block index, explore the random magnetic order. The Parisi matrix Qp ar isi has the 



wellknown form apart from the nonvanishing diagonal elements g; their presence is required by the fact that 
(<r z ) 2 = (n-f — niY 1- The structure of the Parisi-matrix is of course responsible for the rather complicated form 
of the Lagrangian; despite this complication the fermion fields can be eliminated in the standard way, which leads to 
the selfconsistent equations given below. 

It is known since Parisi's work [jjj that an analytical low temperature expansion is hard to obtain even for the standard 
SK-model and its smaller set of selfconsistent parameters. First insight is gained by the one-step RPSB (K = 1). 
The standard three parameter set of the SK-model for K — 1, order parameters qi and q2, and m = m\ ~ T, is 
enlarged in the fermionic space by q — q\ ~ T, where q :— [< <t{t)o{t') >) av represents a spin correlation, which 
remains static unless a fermion hopping mechanism or other noncommuting parts are included in the Hamiltonian. 
For the fermionic Ising spin glass the {K = l)-DoS reads 

p°i E ) = — m ,~ \ 7 -m= — / dy 2 e 92 ttt^2 (3) 

J — oo 1 

with C = cosh{f3H)+Q, where £ = cosh(f3 '^)ea;p(— '(q — qi j) reveals the competition between the particle "pressure" 
exerted by the chemical potential [i and the single-valley susceptibility \ — P(q~Qi) leading to a crossover at \n\ = \\ 
in the T — > 0-limit. The i^-term is a fermionic feature, absent from the standard SK-model. It is closely related to 
the fermion filling; this filling factor behaves discontinuously on the T = 0-axis . For T=0we obtain exactly 
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and a(H) = rn'(T = 0). The replica-symmetric solution displays a magnetic hardgap of width 2E g (H) in the DoS 
and the system remains half-filled at T — within the finite interval given by < ^E g (H). For higher values of the 
chemical potential, hence smaller spin density, phase separation occurs together with a discontinuous transition into a 
full or an empty system pl| . A stable homogeneous saddle-point solution could only be found for the half-filled case. 
Thus the following analysis is restricted to this interval of chemical potentials. Its width is determined selfconsistently 
and seen to decrease to zero as K — > 0. The selfconsistent equations for q,qi,q2 and the Parisi parameter m 0] 
simplify in the T = 0-limit becoming 



q = qi = 1, lirriT^o- 
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where t = q\ — q2- For zero magnetic field one finds q2 — 0.476875, a = HmT-+om'(T = 0) = 1.36104, and \ = 
UiriT^o(q — qi)/T — .239449. The H-dependent solutions shown in fig.l are obtained numerically and then used 
in evaluating eq.([l]) for the density of states. T = 0-results are shown in figs. 2 and 3, while the result at finite low 
temperature of fig. 4 illustrates the presence of plateaus of constant slope, each corresponding to Parisi order parameter 
separations (here: q — qi and q\ — 52)- The number of these plateaus of constant slope increases with the order K 
of Parisi-RPSB. Hence, the time-dependence of the Green's function should characteristically depend on the order 
parameter separations qk — qu-i- If we compare with the replica-symmetric result a reduction of the gapwidth is 
observed. Analytically one finds a gapwidth 



E g (H) =x = lim T ->of3(q - <7i) 



(8) 



which turns into limT->o([3{q — 9(1))) in terms of the Parisi function at K = 00. Only at zero RPSB this susceptibility 
coincides with the equilibrium x- In fact for 1-step RPSB the fermionic Ising spin glass approaches \ — P(q ~ <Zi) + 
(3m(qi — (72) — * -95, the same numerical value as the one for the SK-model. 





FIG. 3. Effect of one step replica symmetry breaking on the fermionic density of states (DoS) for magnetic fields H = 
(curve c: 1RPSB, a: ORPSB) and H/J = 0.6 (d: 1RPSB, b: ORPSB) 

FIG. 4. Low but finite temperature (T = .01J) form of the zero field DoS in 1-step RPSB. 



While the 1-step RPSB provides already a much better approximation than 0-RPSB it is still unstable towards 
higher RPSB. We have therefore extended given equations to arbitrary K. Apart from the K-invariant relation E g (H) = 
X we find a second invariant with respect to K-th RPSB, including K = oo, which is given by 

lim \mv g (H)P«{ E ) = 2* = \ E gi H ) (9) 

The invariant ratio 1/2 is seen analytically by comparing the formulae for the gapwidth and for p(\E\ = E g + 0) 
(both generalized to arbitrary K) in the T — » 0-limit. For each given K, (half) the gapwidth equals the fermionic 
nonequilibrium susceptibility which turns into x = P{Q ~ ?(!)) differing only by exponentially small terms from 
the SK-model result x = — ?(!)) ~ T [Q, where q(l) denotes the Parisi function q[x) at x = 1. Consequently the 
DoS-hardgaps at finite K terminate in a softgap for K — > oo. Note that we did not have to evaluate the T — 0-Parisi 
function q(x) in order to reach this conclusion. Assuming that the relation between gapwidth and x remains valid (at 
least in good approximation) for short-range models, fluctuation effects should harden the gap. This requires further 
analysis. 

A quantity of particular interest in many-body theories is the Ward identity for charge conservation. It shows 
that Parisi symmetry breaking, as observed in the density of states, exists also in vertex functions. The Ward 
identity for the insulating model can be viewed as one of a metallic spin glass at momentum transfer k = 0, ie 
iu> A n (k = 0, e + u), e) — G(e + u>) — Q(e) (fermion momenta suppressed). A„ is the Fourier transformed three point 
function < T T [a\(r)aii (T')hj(O)] > in terms of the fermion operators. Thus the density-part A„ in the Ward identity 
(a current -part Aj emerges for itinerant models and k ^ 0) obeys 

lim^^o u) Umk-^oA% R (k, e + u>, e) = 2iri p(e, {q r - q r -i}) (10) 

and thus shows that the Parisi form of the DoS, depending on all q r — q r —i or on q{x) for K = oo, also enters the 
vertex function. This will also occur in metallic spin glascs, whence diffusive modes and conductivity are expected to 
depend on Parisi symmetry breaking. While we have proved the existence of a spin-glass hardgap at any finite K > 
with 



Sp(E) ~ \E — w(H)\ ,p|>x,lf<oo]; p(E) ~ \E\ X , [K = oo] (11) 

the pseudogap solution at K = oo has a scaling exponent x, which could eventually become different from one and 
remains to be determined. This pseudogap together with x=l would be slightly reminiscent of the exponent found 
for a superconducting glass unitary nonlinear sigma model [^o| . We remark that the pseudogap-solution given for 
the fermionic spin glass model refers precisely to fi = whereas the hardgap-solutions at finite K happened to be 
stable within finite intervals \p\ < x/2 corresponding to half-filling only at T = 0. The regime beyond half-filling, 
identified as the domain of phase separation in the replica-symmetric solution p~l] ] , requires further analysis at T = 
as well as several metallic-, Kondo-type-, superconducting-, and other model extensions. The new method of Fourier 
Transformations in replica space |l2[ is hoped to facilitate further insight into the difficult K = oo-solutions. 
We note that an overlap distribution function for data clustering shown in |l5| ] and interpreted as a pseudo T = 
problem in a classical spin analogy, revealed, apart from the ratio discussed in Eq.(|^), a remarkable similarit y w ith 
the (H = 0)-density of states. It appears interesting to explore pseudo-(T = 0) neural network problems Q as 
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potential classical partners of fermionic spin glasses. 

Summarizing our results we proved i) replica permutation symmetry breaking to be most important in the T = 
quantum field theory of the fermionic spin glass, ii) that low energy excitations are determined by RPSB and hence 
the long-time quantum-dynamical behaviour of the fermion Green's function carries RPSB-fingerprints, which iii) 
affects higher order correlations by means of charge conservation too. 
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to him at the ICTP Trieste (Italy) during the workshop on Statistical physics of frustrated systems, where part of this 
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and David Sherrington. 
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